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Abstract. Johnson noise thermometers (JNT) measure the equilibrium electrical 
noise, proportional to thermodynamic temperature, of a sensing resistor. In the 
correlation method, the same resistor is connected to two amplifiers and a correlation 
of their outputs is performed, in order to reject amplifiers' noise. Such rejection is 
not perfect: the residual correlation gives a systematic error in the JNT reading. 
In order to put an upper limit, or to achieve a correction, for such error, a careful 
electrical modelling of the amplifiers and connections must be performed. Standard 
numerical simulation tools are inadequate for such modelling. In literature, evaluations 
have been performed by the painstaking solving of analytical modelling. We propose 
an evaluation procedure for the JNT error due to residual correlations which blends 
analytical and numerical approaches, with the benefits of both: a rigorous and 
accurate circuit noise modelling, and a fast and flexible evaluation with an user-friendly 
commercial tool. The method is applied to a simple but very effective ultralow-noise 
amplifier employed in a working JNT. 
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1. Introduction 

A Johnson noise thermometer (JNT) determines the thermodynamic temperature T of 
a sensing resistor R by measuring its equilibrium voltage noise, given by the Johnson- 
Nyquist relation [1] S v (f) = Ak-gTR, where S v (f) is the noise voltage spectral density 
function and &b is the Boltzmann constant^. 

Johnson noise thermometry has been employed successfully to measure 
temperatures in a wide range (from below 100 mK to 1800 °C [2]) and to determine 
temperature ratios between fixed points of the International Temperature Scale ITS- 
90 (see [3] for a review). In the framework of a redefinition of the Systeme 
international d'unites (SI) in terms of fundamental constants, the use of JNTs 
performing measurements traceable to the representation of the volt given by the 
Josephson effect has been proposed as a method towards new determinations of k-Q [4-6]. 

Most JNTs employ room-temperature electronics with junction field-effect 
transistor ( JFET) front-end amplifiers. Such amplifiers have an equivalent input voltage 
noise having a spectral density of the same order of magnitude of Johnson noise. 

The use of correlation permits to reject to a large extent amplifiers' voltage noise 
(see [3, par. 6.4] for a review). The resistor noise is amplified using two identical 
amplifiers, and only the correlated part of the amplified signals is considered for the 
measurement. Thanks to advances in mixed-signal electronics and real-time processing, 
digital implementation of correlators is preferred in more recent JNTs. 

The rejection of amplifiers' noise given by the correlation method is not perfect. 
Amplifiers' voltage and current noises are coupled through the sensing resistor and its 
wiring, and give rise to residual correlation terms, indistinguishable from that given by 
Johnson noise, therefore giving a systematic offset in the temperature reading. 

The amount of residual correlation is essentially dependent on the design of the 
input stage of the amplifier and the choice of its components. Some contributions to 
this residual correlation can be determined experimentally, with extreme difficulty, by 
regressions on measurements with different values of R of the same temperature [7]; 
other contributions are beyond direct experimental assessment. 

Therefore, an evaluation must be performed by electrical modelling. We'll see in 
section [2] that the modelling has to take into account not only all relevant voltage and 
current noises generated by different sources, but also the correlations among them. 
Unfortunately, standard packages for numerical circuit simulations, such as SPICE [8], 
do not take into account such correlations, and therefore cannot be employed directly. 

An alternative approach is analytical modelling. Equations modelling all circuit 
elements of the amplifier, and their corresponding noise sources, are written explicitly 
and coupled to circuit network equations [9,10]. The resulting equations can be solved 
analytically. 

The analytical method has been applied to error analysis in JNTs since long [11]; 

| The low-frequency approximation of Johnson- Nyquist relation is reported. The expression is accurate 
to one part in 10 7 at room temperature and frequency below 1 MHz. 
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White and Zimmermann [12] have made a comparative analysis of some simple amplifier 
circuits employed in JNT. The analytical method is complex and heavy, and in practice 
cannot be applied except for the simplest circuits. Any modification in the circuit mesh 
require modifications in the model equations, which must be solved again from scratch. 

In the following, we propose a new method which blends the approach of the 
analytical method with numerical simulation of circuits. The method permits a rigorous 
noise analysis, taking into account noise correlations, with the practical advantages of 
numerical simulations (which an amplifier designer would probably employ in any case). 
The method will be applied to a working amplifier, very simple but of practical interest, 
successfully employed in a working JNT under development [6] . 

The following symbols conventions are used throughout this article: signals are 
represented by lowercase symbols with lowercase subscripts, whereas their Fourier 
transforms are represented by the corresponding uppercase symbols with the same 
lowercase subscripts (we will also consider "Fourier transforms" of random signals: for 
a more thorough discussion see [9, ch. 3]); if x(t) and y(t) are both random signals, 
S x (/) and S y (f) represent their one-sided spectral density functions (/ is the Fourier 
frequency) and S xy (f) represents their cross-spectral density function; for complex 
quantities, a star denotes complex conjugation and the operator Re takes the real part. 

2. Error analysis in unbalanced circuits 

2.1. Error equation 

In unbalanced measurements, a single resistor R generates the noise voltage v(t) which 
is then amplified by two amplifiers with unbalanced inputs and equal gains A v (f). The 
amplifiers' output signals v Q i(t) and v Q 2(t) are then digitized and the cross-spectral 
density function between the two channels is estimated. 



The equivalent circuit for this kind of measurement is shown in figure [TJ For 
simplicity, we will assume that the two channels have identical properties: where 
necessary, the given results can be easily generalized to consider unequal channels. The 
interconnection cables are modelled by the two impedances Z c (f), while Z\(f) represents 
the input impedance of the amplifiers. Voltage sources v a \(t) and v&it) represent the 
equivalent noise voltage of the amplifiers; current sources z a i(£) and -z a2 (t) represent the 
input short-circuit noise current of the amplifiers. We will further assume that noise 
sources belonging to different amplifiers are uncorrelated: with this assumption, only 
the cross spectral densities S Val i al (f) = S Va2 i a2 (f) = S Va i a (f) need to be considered. 

Let T = kftT and 



[Figure 1 about here.] 
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l+2R[l/Z i (f) + l/Z c (f)Y 
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If we assume the product A(f)A v (f) to be known from gain calibration — which is 
typically the case with absolute JNTs [13], then it can be shown that the relative 
measurement error of T due to the amplifiers is given by 

~T = 2T \ RSi » {f) \AjfW J ' (2) 

The above equation is equivalent to (6) of [14] and to (4) of [12]. Generally, the error 
estimated from (T5]) can be frequency- dependent: in such case, an average value in the 
frequency range of interest should be determined. 

2.2. Estimation of the spectral density functions S^(f) and S v ^(f) 

Although it is conceivable to measure <Si a (/) and S v& i & (f), in practice such operation is 
difficult and cannot be done with sufficient accuracy over the whole frequency range of 
interest. Therefore, some kind of analysis should be employed to evaluate S^(f) and 
S Va i a (f). In [12], the analysis has been done with an analytical method; however, this 
approach is cumbersome and error-prone, and could only be employed with the simplest 
circuits. To circumvent these problems, we propose a method based on numerical 
simulations which is applicable to more complex circuits. 

Consider one of the two amplifiers and let Xj(t), which can be either a voltage 
or a current, be the jth noise source internal to the amplifier. Assuming the input 
short-circuited, the equivalent input noise voltage v & {t) and the input short-circuit noise 
current i a (t) of the amplifier can be written as a superposition of the noises generated 
by each source Xj(t), i.e. 

W) = = £<&(/)*>(/) (3a) 
/.(/) = 5>i(/)x>(/) (3b) 



j 



The network functions 



Q. {f)= Wl = 1 W) (4a) 
W) - |§ (40) 



can be approximately estimated by means of an analog circuit simulation prograrnj. In 
particular, such programs can perform small-signal analysis to determine directly the 
ratios V (f)/Xj(f), I a (f)/Xj(f) and A v (f) and the input impedance Z^f). 

Using ( |3fflj) -(l3fej). it is not difficult to show that the searched spectral densities can 
be written as (S Va (f) is given for completeness, although it does not enter explicitly 
in ©) 

SJ/) = £ty(/)«(/)S*«*(/) (5a) 

§ In this work, we have used LTspice IV® from Linear Technology Corporation [15]. 



Systematic errors in JNTs 



5 



SiM)-J2 P ^ P ^f) S ^M) (56) 

W/) = J2Qj{f)Pk(f)S* i * k (J) (5c) 

where S XjXh (f) is the cross-spectral density function between the noise sources Xj(t) and 
Xk(t), with the convention that, if j = k, S XjX . (/) represents the spectral density function 
of Xj(t). The assignment of the functions S XXk (f) must be based on the known noise 
models of the devices composing the specific amplifier under examination. Actually, 
most of the S XjXk (f) with j ^ k can be considered identically zero. 

Although not easy to evaluate, the uncertainties of the estimations given by 
equations (!5aj) - (l5cj) are related to the uncertainties of the various parameters on which 
the functions Pj{f), Qj(f) and S XjXk (f) depend, and would be the same obtainable with 
the analytic method. 



3. Error analysis in balanced circuits 

Typically, in order to improve the rejection to external disturbances, balanced 
measurements are preferred over unbalanced ones. The equivalent circuit for a balanced 
measurement is shown in figure EJ two instrumentation amplifiers with differential gain 
A d (f) measure the differential noise voltage generated by two sensing resistors working 
at the same thermodynamic temperature T. Again, the amplifiers' output signals are 
digitized and the cross-spectral density function between the two channels is estimated. 
The symbols in figure [2] have the following meaning: v ip (t) and v in (t) are the noise 
voltages generated by the sensing resistors R; Z cc (f) models the impedance of the 
interconnection cables; Z id (f) and Z ic (f) represent respectively the differential-mode 
and common-mode input impedances of the amplifiers [22]; f a i(t) and v^if) represent 
the equivalent noise voltage of the amplifiers; i a ip(*)) i&init), ia.2 P {t) and i a 2n(t) represent 
the input short-circuit noise currents of the amplifiers. 

[Figure 2 about here.] 



Let 



A(f) _ E dl (f) _ E d2 (f) 

1 (66) 



l + 2R[2/Z iA (f) + l/Z cc (f)} 
then, in this case, with the same assumptions made in Sj2j the relative measurement 
error of T = k^T due to the amplifiers is given by 



AT 1 



T 4T 



2R - Re S^(/)] 



Re i A (f) s^jn-s^jf) 
\MfW 



(7) 
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where S{ & (f) is the spectral density function of iai p (t), i&in(t), i a 2p(t) and i a 2 P (t), 
Kaplan (/) * s the cross-spectral density function between i a i p (t) and i a i n (t) (or between 
ia,2p(t) and i a 2 a (t)), S v& i ap (f) is the cross-spectral density function between v a i(t) and 
i a ip(t) ( or between v^it) and i a 2 P (t)), ^^(J) is the cross-spectral density function 
between t» a i(i) and i a i n {t) (or between v a 2(f) and ^ a 2n(^))- 

The procedure outlined in ^2] can be generalized to yield the above defined spectral 
densities: 

SJJ) = £Qi(/)%(/)<W/) (8a) 

3>k 

SiM) = J2 p M p vk(f)s^ k (f) = Y, p ^) p ^f)s XjXk U) (86) 

^apian(/) = £ ^pi(J) ^ (/) S XjXk (/) (8 C ) 

S^(f) = E^^) P n fc (/)^(/) (8e) 



with 



/>„(/) = |$ (96) 

Also in this case, an analog circuit simulation program can be used to estimate the 
network functions Qj(f), Ppj(f) and P n j(f). 



4. Example of application 

The INRIM JNT [6] employs two amplifiers with the circuit of figure [31 Despite its 
simplicity, it has a voltage noise below InV/VThz at audio frequencies, which is much 
lower than the best JFET integrated circuit amplifiers, and comparable to ultralow- 
noise discrete amplifiers (see e.g. [16-21]) which are however more complex. The value 
of the sensing resistor is 1 kfl and the interconnection cables have a capacitance of about 
80 pF. 

[Figure 3 about here.] 
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4-1. Description 

Looking at figure [3] we see that the JFET Q, directly connected to the input voltage Vi, 
works in a common-source configuration with a gate-source bias voltagqjj] Vgs ~ V. 

The operational amplifier OA works as a transresistance amplifier with gain Rf, and 
biases Q; the transresistance configuration eliminates the Miller effect [22, ch. 7], thus 
enhancing the bandwidth. The drain-source voltage Vds is set by OA at the voltage 
Vb ~ 9 V of the polarization battery, which has an extremely low noise [23,24]. 

OA works in a single-supply configuration; its output v Q is ac-coupled through 
capacitor C, and can be further amplified by additional stages if necessary. 

The overall low-frequency gain A v0 = g m Rf of the amplifier depends on Rf and 
on the transconductance g m of the JFET; the bandwidth of the amplifier depends on 
the JFET capacitances and on the gain-bandwidth product of OA; such parameters 
can have significant deviations from sample to sample and have a strong temperature 
dependence. In JNT applications, where periodic calibration [3] of A v (f) is performed 
with a reference signal, this does not constitute a problem. 

Typical supply voltage is Vg ~ 24 V from an unregulated battery; the quiescent 
current is about 15 mA. 

4-2. Construction and performance 

A prototype has been assembled with a 2SK170 JFET (Toshiba), an OP27 (various 
suppliers) for OA and Rf —1MI. The gate current I g is 2pA-h3pA (measured with a 
Keithley mod. 6430 current meter), which gives a current shot noise less than 1 fA/ vTfz. 

The transfer function of the amplifier is shown in figure HI It has been measured 
with a network analyzer (Agilent Tech. mod. 4395A), injecting the signal with a resistive 
divider (50 Sl-0.5 O). A 3 dB-bandwidth of about 4 MHz can be estimated, with a gain 
flatness better than 1 dB up to 1 MHz. 

[Figure 4 about here.] 

The square root of S Va (f) is shown in figure EJ It has been measured with a two- 
channel signal analyzer (Agilent Tech. mod. 35670A) by connecting v Q to both channels 
and performing a cross-correlation measurement in order to reject the analyzer noise. 
The noise floor is about 0.8nV/VTIz, corresponding to the Johnson noise of a ^38 Vl 
resistor at 300 K. We do not consider the noise spectrum for frequencies below 1 kHz, 
and hence we neglect the effect of low-frequency noises such as generation-recombination 
and flicker noises, because the operating frequencies of JNTs are typically well above 
1kHz. 

[Figure 5 about here.] 

|| A drawback of such configuration is the limited dynamic range, a few tens of mV at the input, before 
distortion occurs; however, in a JNT the integrated signal is of a few pV. 
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4-3. Results 

[Figure 6 about here.] 

Figure [6] shows the small-signal equivalent circuit of the amplifiers with the five noise 
sources considered in the error analysis: v\{t) and i<i (t) represent respectively the 
equivalent voltage and current input noise generators of the JFET Q; i 3 (t) represents 
the thermal noise of Rf, i±{t) and v$(t) represent respectively the equivalent current and 
voltage input noise generators of the operational amplifier OA. V\(t) and i 2 (t) have been 
assumed with spectral density functions [9, ch. 5] (C gs is the gate-source capacitance of 
the JFET) 

3.9m 

^ 2 (/)^2g/ g + (27r/Cgs)2 fc B T (106) 
and cross-spectral density function [9, ch. 5] 

^(/)«0.4j [S Vl (f)S l2 (f)} 1/2 ■ (10c) 

is(t) has been assumed with spectral density function Si 3 (f) = Ak^T/R^ while the 
spectral densities of i&{t) and v^(t) have been inferred from the operational amplifier's 
data sheet. Apart from Vi(t) and ^(t), all other generators have been considered 
uncorrelated. 

[Figure 7 about here.] 

The expected temperature measurement error resulting from the simulation is 
shown in figure as a function of frequency. For low frequencies the error tends 
asymptotically to a constant: here, the main error source is the gate current shot noise of 
the JFET. At high frequencies the error tends to increase as f 2 : this behaviour depends 
mainly on: i) S i2 (f) which increases as f 2 , ii) t>i(t) which contributes to ijt) through 
the input impedance of the amplifier, and iii) v^{t) which contributes to i a (t) through 
the gate-source capacitance of the JFET. 

From figure [7J the error can be approximated by 
AT 

=-«1.2x 1(T 7 + 3.04 x 10- 14 (//Hz) 2 (11) 

and the expected average error over the present operating frequency range of the INRIM 
JNT (3kHz^-7kHz) is about 9.2 x 10~ 7 . 



5. Conclusions 



We have described a method for the evaluation of the systematic error of the temperature 
reading of a correlation JNT caused by residual terms in the correlation spectrum of the 
amplified Johnson noise. Although a general expression for such error can be given 
(equation (j7])), its terms are dependent on circuit components and topology. The 
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method consists in the identification of noise sources of the components employed in 
the amplifiers' network; in the use of a commercial circuit simulation software, with its 
user-friendly interface, to numerically estimate a number of transfer functions for such 
noise signals; and in a numerical calculation (equations (|8aj) - (l9cj) ) easy to implement 
on any platform. Versions of the model for both single-ended and differential input 
amplifiers are given. 

The paper shows the application of the method to a simple but effective single- 
ended open- loop amplifier of practical interest. However, the method can be easily 
applied to much more complex amplifier networks for which an analytical approach 
would be unfeasible. 

Future developments will be devoted to the application of the method to other 
JNT amplifiers, either already working or under development, having both open- and 
closed-loop topologies. 
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Figure 1. Equivalent circuit used in the analysis of unbalanced measurements. 
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Figure 2. Equivalent circuit used in the analysis of balanced measurements: Z- 1( \ 
represents the differential-mode input impedance of the amplifiers; Z lc represents the 
common-mode input impedance; and || (— 2Z- 1C ) = —2Z XC Z^I{Z^ — 2Z lc ) (see [22] 
for details). 
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Figure 3. The amplifier schematics. See text for an explanation of the symbols. 



FIGURES 



15 



40 I rTTT| ,— , 



50 




"| 5 l 1 1 1 1 1 i i i 

1kHz 10kHz 100kHz 



50 -0 



0) 

w 

CD 



100^ 



* -150 



-200 



1MHz 10MHz 
Frequency 



Figure 4. The transfer function of the amplifier, in magnitude and phase 
representation. 



FIGURES 




Figure 5. Equivalent input voltage noise of the amplifier. 
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Figure 6. Small-signal equivalent circuit of the amplifier of figure[3]with noise sources. 
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Figure 7. Expected temperature measurement error as a function of frequency for 
the amplifier of figure [3] (more details in text). 



